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Abstract. We present observations of the stretching 
mode of ^■^C02 ice along 13 lines of sight in the Galaxy, 
using the Short Wavelength Spectrometer on board of the 
Infrared Space Observatory (ISO-SWS). Remarkable vari- 
ations are seen in the absorption band profile in the differ- 
ent lines of sight. The main feature is attributed to ^^C02 
mixed with polar molecules such as H2O, and CH3OH. 
The high-mass protostars GL 2136, GL 2591, S 140 : IRSl, 
and W 3 : IRS5 show an additional narrow substructure 
at 2282 cm~^ (4.382 /im), which we attribute to a po- 
lar, CHaOH-containing CO2 ice, that experienced heat- 
ing. This heating effect is sustained by a good correlation 
of the strength of the substructure with dust and CO gas 
temperatures along the line of sight, and anti-correlation 
with ice abundances. Thus, our main conclusion is that 
interstellar CO2 ices around luminous protostars are sub- 
jected to, and altered by, thermal processing and that it 
may reflect the evolutionary stage of the nearby protostar. 
In contrast, the ices around low mass protostars and in a 
quiescent cloud in our sample do not show signs of thermal 
processing. 

Furthermore, we determine for the first time the Galac- 
tic ^^C/^^C ratio from the solid state as a function of 
Galacto-centric radius. The ^^C02/^'^C02 ratio for the 
local ISM (69±15), as well as the dependence on Galacto- 
centric radius, are in good agreement with gas phase 
(C^^O, H2CO) studies. For the few individual objects for 
which gas phase values are available, the ^^C/^^C ratios 
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derived from CO2 tend to be higher compared to CO stud- 
ies (albeit with ^ 2.5 a significance only). We discuss the 
implications of this possible difference for the chemical 
origin of interstellar CO2. 

Key words: ISM: dust, extinction - molecules - abun- 
dances - Infrared: ISM:lines and bands - Galaxy: abun- 
dances - Stars: formation 



1. Introduction 

Space based observations have revealed that the CO2 
molecule is an important constituent of quiescent and 
star f orming molecular cl ouds (d'Hendecourt & de M uizon 

d'Hendecourt et al 



1989: de Graauw et al. 1996 



Giirtler et al. 



1998; Gerakines et al 



1996 



199q; Str azzul la et al. |1998| ; Whittet et al. 

1999| ). It is primarily present in 



the solid state (Van Dishoeck et al. 1996), with an abun- 
dance of typically 20% of water ice, consuming ~3% of the 
carbon budget. The main focus of research has been on 
the high quality observations of the solid ^^C02 stretch- 
ing and bending modes at 4.27 and 15.2 /im, which are 
prominently present in the spectra obtained with ISO- 
SWS (Gerakines et al. |1999| ; Boogert |1999| ). 

The detection of the ^^C02 isotope, with a two or- 
ders of magnitude lower abundance, has been reported as 
well (de Graauw et al. 1996; d'Hendecourt et al. 1996). Its 
stretching mode frequency is well separated from that of 
^^C02 (4.38 fim vs. 4.27 /im). The analysis of this weak 
band is very attractive and has some specific advantages 
over ^^C02 studies. It is an independent and very sensi- 
tive probe of the ice mantle composition. "'^^002 is always 
a trace constituent. In contrast to ^^C02, interactions of 
^^C02 molecules with themselves, and consequent effects 
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on the band profile, are unimportant. For the same rea- 
son, the polarizability of the ice at the ^^C02 resonance 
frequency can be neglected, and the absorption profile is 
insensitive to the shape and thickness of the ice mantle, 
and composition of the grain core, contrary to the ^^C02 
bands. Thus the laboratory spectrum can be compared 
directly to the interstellar spectrum, and uncertainties re- 
sulting from corrections induced by the grain shape or 
uncertain optical constants are avoided (Ehrenfreund et 
al. p^97| ). 

The main motivation for our study is to determine the 
physical and chemical history of interstellar CO2, and in- 
terstellar ices in general. Once interstellar ices are formed, 
by a combination of direct accretion from the gas phase 
and chemical reactions on grain surfaces, they can be ex- 
posed to a variety of violent processes. Among these are 
cosmic ray bombardment, ultraviolet irradiation, heating 
by visible and infrared photons, and disruption in shocks. 
What is the importance of these processes in various en- 
vironments, such as protostars of low and high mass, and 
quiescent clouds? Can the ices around low mass proto- 
stars survive the various destruction mechanisms, and be 
included in comets? If so, how much do its composition 
and structure still resemble the interstellar ices? In or- 
der to answer these questions, we analyze the absorption 
band profile of the stretching mode of interstellar ^^C02. 
We make use of a large database of spectra of CO2 ices, 
with a range of compositions and temperatures, obtained 
in the Leiden Observatory Laboratory (Ehrenfreund et al. 



1996 ; 1997 ; 1999 ). In a separate study, the analysis of the 
absor ption bands of ^^C02 will be presented (Gerakines 



et al. |1999D . 

Another important motivation for this study is to de- 
termine the ^^C02 abundance, and derive the Galactic 
i2(-;^i3(^ abundance ratio. It is the first time that the 
i2(-;^i3(-; j-g^^Q jg determined from the solid state. An im- 
portant advantage over gas phase studies is that the col- 
umn density can be straightforwardly derived from the 
observed integrated optical depth and the intrinsic band 
strength determined in the laboratory. In contrast, for gas 
phase species excitation, and radiative transfer models are 
required. 

Previous (gas phase C^'^O, II2CO) studies have shown 
that the ^^C/^^C ratio increases with Galacto-centric 
radius, with ^^C/^'^C=25 in the Galactic Center, and 
i2(-;^i3(^ _ ^Yie local interstellar medium (see Wil- 



son & Rood 1994| for an overview) . Recent observations of 
atomic C and C+ y ield ^^C/^'^C=60 toward the Orion Bar 
(Keene et al. 1998|) . Determination of the ^^C/^^C ratio is 



an important input for evolutionary models of our Galaxy, 
since ^^C is produced by Helium burning by massive stars, 
which can be converted to ^^C in the CNO cycle of low- 
and intermediate-mass stars at later times. 

Furthermore, comparison of the ^'^C/^^C ratios derived 
from various species, will allow to determine the impor- 
tance of chemical fractionation ("'^'^C preferentially incor- 



porated in CO) and isotope-selective destruction (^^CO 
preferentially destroyed). Models of photo-dissociation 
regions (PDRs), including chemical fractionation and 
isotope-selective destruction, show that the ^^C/^^C ra- 
tio for gaseous C^^O can decrease by '^50%. Recent ob- 
servations of C and C^, however, indicate that chemical 
fractionation is not an important effect, or is compensated 
for by isotope-selective photo dissociation, in these PDRs 



(Keene et al. 1998). We do not expect that these effects 
play an important role for the CO in our lines of sight, 
since they mainly trace dense molecular cloud material 
with low atomic C abundances. However, species that are 
formed from atomic C or C^ rather than CO, ma y hav e 
very different ^^C/^^C ratios (e.g. H2CO; Tielens |l997| ). 
Thus, the determination of the solid ^^C02/^^C02 ratio is 
a potentially powerful tool to trace the chemical history of 
interstellar CO2; does it originate from CO, as is generally 
assumed, or perhaps from C^"*"^? 

In Sect. 2 we present our ISO-SWS observations and 
data reduction techniques. The laboratory results of the 
^^C02 stretching mode are summarized in Sect. 3. The 
analysis of the interstellar ^^^002 band profile, and deriva- 
tion of the column densities is presented in Sect. 4. We 
discuss the main results in Sect. 5. First, we correlate the 
depth of the detected narrow 2282 cm~^ absorption fea- 
ture with known physical parameters along the observed 
lines of sight. Then we compare the ^^C/^'^C ratios, de- 
rived from solid CO2 and gas phase molecules, and discuss 
the importance of chemical fractionation, and the reaction 
pathway to form CO2. The summary and conclusions are 
given in Sect. 6. 

2. Observations 

The "CO2 spectra were obtained with ISO-SWS (de 
Graauw et al. 



1996 



1996| ). Most spectra 



Kessler et al. 

were observed in the high resolution grating mode (i? 
A/AA=1500; "SWS06" mode). Only the source GL 989 
was observed in the fast scanning mode ("SWSOl speed 
3") at an effective resolving power R ~400. 

The spectra were reduced with version 6 of the SWS 
pipeline during January- April 1998 at SRON Groningen, 
using the latest calibration files available at that time. For 
some low flux spectra we applied a "p ulse sh ape" correc- 
tion on ERD level (de Graauw et al. 1996) to linearize 
the slopes, resulting in ^^10% better signal-to-noise. All 
detector scans were checked for excessive noise levels, de- 
viating flux levels and continuum slopes, and dark current 
jumps. Bad scans were taken out off the data, and then 
the scans were "flat-fielded" to the median spectrum us- 
ing a first order polynomial. Cosmic ray hits were removed 
by clipping all points deviating more than 3 a from the 
median, and subsequently the scans were averaged and re- 
binned per scan direction to a resolving power of i? = 1500 
(R = 400 for GL 989). The effective "signal-to-noise" 
values, as determined from the systematic differences be- 
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Table 1. Observed lines of sight 



Object 


RA (1950.0) 


Dec (1950.0) 


Revolution" 


km s^^ 


S/N" 
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497 


— 4'H 




Plins Ifi 






uou 


-1-1 fi 


vo 








UU 1 


-t-97 


UiJ 


GL 989 


06:38:25.3 


+09:32:29 


716 


d 


100 


Elias 29 


16:24:07.7 


-24:30:40 


452 


-7 


100 


GC 3 


17:43:04.4 


-28:48:27 


327 


e 


140 


W 33A 


18:11:44.2 


-17:52:59 


467 


+20 


65 


GL 2136 


18:19:36.6 


-13:31:40 


872 


+13 


100 


R CrA : IRS 2 


18:58:19.0 


-37:02:50 


495 


0.0 


50 


HH 100 


18:58:28.2 


-37:02:29 


704 


0.0 


65 


GL 2591 


20:27:35.8 


+40:01:14 


357 


-25 


200 


S 140 : IRSl 


22:17:41.1 


+63:03:42 


263 


-20 


200 


NGC 7538 : IRS9 


23:11:52.8 


+61:10:59 


433 


-70 


100 



^ISO revolution number 

''Heliocentric velocity from literature 

'^Effective signal-to- noise of the '^C02 spectrum 

'^No radial velocity correction applied; spectrum has low resolution 

°No radial velocity found in literature (see text) 



tween the SWS up and down scans, are given in Table 
We note that due to the weakness of the interstellar ^^C02 
band, ~ 10% of the continuum, our analysis (column den- 
sities, band profiles) is not significantly affected by de- 
tector memory effects, known to exist at this wavelength. 
Also, the detector response function is very smooth over 
the narrow ^^C02 band, and thus will not induce spuri- 
ous features due to small wavelength calibration errors or 
dark current uncertainties. 

An accurate wavelength calibration is important for 
our study. The standard SWS wavelength calibration was 
applied (Valentijn et al. 1996). The satellite velocity was 
converted to the heliocentric reference frame, using stan- 
dard ISO-SWS pipeline routines. Then we corrected for 
the source velocity in the heliocentric frame, using pub- 
lished radial velocities for each source (Table |l|). Recent 
millimeter observations of C^^O emission lines show veloc- 
ity shifts that are in excellent agreement with the values 
that we use (Van der Tak et al., in prep.). To check the 
accuracy of the wavelength calibration and velocity cor- 
rections, we compared the position of ro- vibrational ^^CO 
gas phase lines (^^CO lines are too weak), when present 
in the same spectrum, with peak positions from the HI- 
TRAN database (Rothman et al. |1992D . For S 140 : IRSl, 
and W 3 : IRS5 the lines are blue-shifted by 10% of a 
resolution element (~0.15 cm~^; ^20 km s~^) compared 
to the HITRAN wavelengths, for GL 2591 by 20%, and 
for GL 2136 by 30%. Part of this blue-shift may be real. 
Mitchell et al. ( |1991[ ) find that, whereas the velocity shift 
of ^'^CO lines is in good agreement with millimeter stud- 
ies, the ^^CO absorption lines show significant blue shifted 



shoulders and sub-peaks, originating from the outfiow as- 
sociated with the protostar. At the much lower resolution 
of our observations ('-^^160 km s^^), this may result in a 
blue shift of the ^^CO absorption line. However, the main 
absorption is present at the cloud velocity, and at our res- 
olution the peak would shift by not more than —15 km s~^ 
(GL 2591), -10 km s^^ (W 3 : IRS5) or much less than 
that for the other sources. Thus, we attribute the CO line 
shifts in the ISO-SWS spectra primarily to uncertainties 
in the wavelength calibration and small pointing errors. 
In particular, GL 2136 was observed at the end of the ISO 
mission for which no updated wavelength calibration files 
are available yet. We applied the small shift derived from 
the CO lines to further improve the wavelength calibra- 
tion of the "'^'^002 band. For one object, GC 3, we found 
no radial velocity observation in the literature, and deter- 
mined it from the CO lines to be WhcUo = -90 km s-i. 
The fully reduced spectra are given in Fig. [T]. 

3. Laboratory CO2 studies 

The absorption band profiles of the stretching and bending 
modes of ^^C02 and 1^002 in apolar ices (CO, N2, O2) 
were extensively studied in the laboratory by Ehrenfreund 
et al. (1996; 1997). Additional experiments of CO2 mixed 
with polar molecules (H2O, CH3OH) were performed in 
the Leiden Observatory Laboratory as well (Ehrenfreund 
et al., 1998; 199E). Figure ^ summarizes the peak shifts 
and broadenings of the stretching mode of ^^^002 in these 
ices at a temperature T =10 K. 

For a pure CO2 ice, i.e. 12002:^3002=90:1, the ab- 
sorption band of ^^002 is centered on 2283.0 cm"^ and is 
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Fig. 1. The ISO-SWS spectra of ^'^C02 on optical depth scale toward massive protostars (left two panels). The right 
panel shows low-mass protostars, a background star of the Taurus Molecular Cloud (Elias 16) and a Galactic Center 
source (GC 3). The spectra were normalized to the peak optical depth by multiplication with the factors given in 
square brackets. The error bars were determined from the difference between the ISO-SWS up and down scans, and 
thus include systematic errors. The vertical dashed lines are given to facilitate comparison of the spectra in the different 
panels. 



very narrow (FWHM~2.3 cm-i). A large broadening and 
shift to lower frequencies is observed when II2O is added. 
An even larger peak shift, up to 2274 cm~^, but a small 
narrowing, occurs in C02:CH30H ices. This difference in 
spectroscopic behavior of Il20:C02 and CH30H:C02 ices 
is also particularly evident in the ^^C02 bending mode 
(Ehrenfreund et al. 1998 ; 1999 ). It is ascribed to the for- 
mation of stable CII3OH.CO2 complexes. Thus, at low 
temperatures, H2O and CHaOH-rich CO2 ices have very 
distinct ^"^002 peak positions. Mixtures of CO2 with both 
H2O and CII3OII lie in between these extremes, as in- 
dicated by the "weak" (H2O:CH3OH:CO2=7:0.6:l) and 
"strong" (H2O: CH3OH: 002=1.7:0.6:1) interstellar mix- 
tures in Fig. 1^. Apolar molecules, such as CO and O2, 
induce much smaller broadenings and peak shifts to the 
^■^002 ice band. Thus, the laboratory simulations show 
that the ^^C02 stretching mode is a sensitive probe to 
discriminate between interstellar polar and apolar ices. 

The peak position and width depend strongly on the 
temperature of the ice. Upon warm-up, the peak posi- 
tion of ^^002 in H2O-, and CIl3 0II-rich ices (with less 
than ~20% of CO2) shifts significantly to lower frequen- 
cies, and at the same time becomes narrower (Fig. |^). 
For example, in the mixture Il2O:CO2=10:l, the ^"^002 
band shifts by 2.5 cm~^, and the FWHM decreases with 
2.5 cm~^, when heated from T =10 to 140 K in the lab- 
oratory. Unfortunately, heating of the polar, C02-poor 



ice has the same effect on the band profile as increasing 
the H2O/CO2 and CH3OH/CO2 mixing ratios (Figs. |a 
and ^). Thus, the ^'^002 stretching mode is less suited 
to determine an accurate II2O/CII3OII/CO2 mixing ratio, 
or the precise temperature of these particular ices. How- 
ever, for H20:C02 and CH30H:C02 ices with CO2 con- 
centrations between 20-90%, heating induces a very spe- 
cific spectroscopic signature. The ^^002 profile becomes 
asymmetric at a laboratory temperature of T '^50 K, and 
^105 K respectively (Figs. || and ^; see Sect. 5.1 for the 
corresponding temperatures in interstellar space). A sec- 
ond peak develops at the blue side of the band, which is 
reflected in Fig. ^ as a "broadening" of the overall pro- 
file. At higher temperatures, this new peak dominates the 
^^C02 spectrum. Its peak position of ~2282-2283 cm~^, 
and narrow width (~3 cm~^) are close to the values for 
a pure CO2 ice, although, in particular for CH3OH ices, 
this correspondence is not exact. This spectroscopic be- 
havior is attributed to the destruction of bonds between 
the polar molecules and CO2, and the formation of new, 
stronger, hydrogen bonds between the polar molecules 
(Ehrenfreund et al. 199£). The CO2 molecules now in- 



teract primarily with themselves. This re-arrangement of 
bonds ("segregation") is also particularly well traced in 
the ^^002 bending mode, which shows, after heating, the 
typical double peak ed structure of a pure CO2 ice (Ehren- 
freund et al., 1999). In accordance with heating experi- 
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Fig. 2. a and b. Laboratory spectroscopy of the stretch- 
ing mode of ^'^C02 ice at T =10 K (a), and at higher 
temperatures (b). The average Gaussian parameters of 
the interstellar spectra are given by crosses. In panel a, 
the peak position versus width variations are shown for 
CO2 in series of increasing H2O, CH3OH, and CO abun- 
dances (connected with a line) . "W" and "S" indicate the 
weak and strong interstellar mixtures (see text). Panel b 
clearly illustrates the large effect of temperature on width 
and peak position of the ^'^C02 stretching mode, due to 
the segregation of CO2. 



ments of H20:CH30H ices (Blake et al. 1991), it seems 
likely that also a spatial segregation exists between the 
CO2 molecules and polar species. Crystalline phases of 
H2O ice enclosing 0.1 /xm size pockets of CH3OH were ob- 
served in these heated H20:CH30H ices. Such microscopic 
observations are essential to determine the structure of 
heated polar CO2 ices as well. 

Finally, for apolar ices containing CO, O2, and N2, 
the ^^C02 band shifts to higher frequencies, and becomes 
narrower upon warm-up. This can be attributed to evap- 
oration of the host molecule. These laboratory results will 
be used to derive the composition of interstellar ice, i.e. 
polar versus apolar, and possibly its temperature history. 



4. Results 
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Fig. 3. Effect of heating on the ^'^C02 stretching mode in 
the ice mixture H2O:CH3OH:CO2 = l:0.92:l. This shows 
clearly the segregation of CO2 at high temperatures. The 
peak at ~2282 cm"'^ is close to the peak position for a pure 
CO2 ice (Fig. ^). We note that the indicated temperature 
was measured in the laboratory. At the much longer time 
scales in interstellar space, the segregation temperature is 
much lower (~78 K; Sect. 5.1). 



4.I. Contamination by gas phase lines 

The spectral region of the ^"^002 ice band is contaminated 
by narrow absorption lines of gas phase CO in the lines 
of sight of GL 2136, GL 2591, and possibly S 140 : IRSl 
(Fig. They originate from CO in the ground vibra- 
tional state (v=0-l), at high rotational levels, up to per- 
haps J=45. The presence of hot CO gas in these lines of 
sig ht wa s shown in other studies as well (e.g. Mitchell et 
al. |1990| ; Table |). The main component of the ^^C02 
ice band is much broader than the unresolved gas phase 
CO lines, and our conclusions are not influenced by this 
contamination. However, the narrow line at 2282 cm~^ de- 
tected in these sources (Sect. 4.2) is at the same frequency 
as the CO J=46^47 R-branch line. It is unlikely that this 
line can be attributed to gaseous CO, given its large depth 
relative to neighboring R-branch lines, in particular for 
S 140 : IRSl. Furthermore, we checked for the presence of 

features from vibrationally excited ^■^CO. The band head 

-1 



of the v=l-2 transition at ~2298 cm"^ (Goorvitch |1994D 
may be present toward GL 2136, and GL 2591 (Fig. ^. 
This feature is well separated from the ^^^002 ice band. 
There is no indication for absorption at the frequency of 
the v=2-3 band head (~2268 cm~^). Given the weakness 
of these ^^CO band heads, the contribution of the v=0-l 
band head of ^^CO (~2276 cm"!) to the "CO2 ice band 
can probably be neglected. 
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One object in our sample, Elias 16, is an evolved 
star located behind the Taurus molecular cloud. It is a 
K 1 giant (Elias 197§| ), and the 2-5 yum spectral region 
shows many narrow absorption lines, most likely asso- 
ciated with the photosphere of this object (Whittet et 
al. 199S). In the spectral region for this study (~2260- 



4.36 



X (wm) 
4.40 



4.44 



2300 cm~^), we detect narrow absorption lines at 2268, 
and 2297 cm~^ (Fig. p. These lines are well separated 
from the region where CO2 ice absorption is expected, 
and are attributed to the v=2-3 and v=l-2 band heads 
of ^^CO (Goorvitch 1994). However, another narrow line 
at 2276 cm-\ probably the v=0-l band head of "CO, 
might be blended with the ^^C02 feature. We compared 
the spectrum of Elias 16 with the K 1.5 giant a Bootis 
(Fig. H). The spectral type of these stars is very similar, 
and indeed a Bootis does show absorption lines at 2268, 
2276, and 2297 cm~^. However, the lines in a Bootis have 
rather prominent wings on the red side, which are not 
recognized in EHas 16. Perhaps the spectral type or lu- 
minosity class of Elias 16 deviates from a K 1 giant. To 
correct for photospheric absorption, we subtracted the op- 
tical depth spectrum of a Bootis from Elias 16 (Fig. 
The narrow lines at 2268, and 2297 cm~^ are nicely re- 
moved, giving confidence in the correctness of this method. 
A significant absorption remains at ~2278 cm~^, which we 
attribute to solid ^^C02 in the Taurus molecular cloud. In 
our further analysis, we used this corrected spectrum of 
Ehas 16. 

4-2. Interstellar solid '^'^002 band profile analysis 

The ^'^C02 band profile varies remarkably between the 
observed lines of sight (Fig. p. Two independent com- 
ponents are present. The high-mass protostars GL 2136, 
S 140 : IRSl, W 3 : IRS5, and GL 2591 show both a 
"broad" absorption extending over 2270-2287 cm~^, and 
a very narrow feature at ~2282 cm~^. Two lines of sight do 
not show the narrow component (Elias 16, GC 3), while in 
others it may be blended with the blue edge of the broader 
feature. In NGC 2024 : IRS2 and GC 3 the profile analysis 
is limited by a low depth and poor signal-to-noise, while 
in GL 989 the resolving power is too low to resolve the 
narrow feature, if present. 

We followed two different approaches to analyze the 
interstellar ^^C02 band profile. First, analogous to the 
solid CO band (e.g. Tielens et al. 1991), the narrow and 
broad component may originate from different ices along 
the line of sight. The ice composition and structure de- 
pend on the chemical (e.g. the H/CO ratio) and physical 
(e.g. temperature) history, which may have varied along 
the line of sight. Therefore, we decomposed the absorp- 
tion band by fitting Gaussians to the observed broad and 
narrow components, and we compared the peak positions 
and widths with the laboratory results (Figs. ||, and 
In our second approach, we investigate whether the ^"^002 
ice band profile can be explained by a single ice. Here, we 
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Fig. 4. Comparison of gas phase CO frequencies with 
sources showing fine structure. The thick vertical lines in- 
dicate the frequency of the R-branch lines of gaseous CO, 
originating from r otational levels J = 29 (2241.7 cm~i) 
to J = 53 (2295.1 cm~^) in the lowest vibrational state 
(v=0-l). The * symbols on top indicate the expected fre- 
quencies of the CO band heads (see text). We attribute 
the sharp feature at 2282 cm~^ and the underlying broader 
component to interstellar solid ^^C02. 

determine the observed peak position and width (not fit- 
ting Gaussians) and again compare these with the labora- 
tory results in Figs. ||, and ^. Both methods provide good 
insight in the characteristics of the observed band profile, 
and the variations between the various lines of sight. 

4.2.1. Two component ices 

For five objects in our sample, good fits are obtained 
by fitting two Gaussians: NGC 7538 : IRS9, GL 2136, 
S 140 : IRSl, GL 2591, and W 3 : IRS5 (Table |). For the 
other lines of sight, no significant improvement is achieved, 
when using two instead of one Gaussian. Nevertheless, 
for some sources the narrow 2282 cm~^ feature might be 
blended with the blue edge of the broader component, and 
we do perform a fit with two Gaussians as well. 

For the broad component we find a Gaussian FWHM 
of typically 6-9 cm~^, with a peak position varying be- 
tween 2276-2280 cm^-'^ for the different objects. Thus, this 
feature can be ascribed to absorption by ^'^C02 ice mixed 
with molecules with a large dipole moment, i.e. H2O and 
CH3OH (Figs. I and ||b; Table 2). To derive the chemical 
history of interstellar CO2 it is important to determine the 
relative amounts of H2O, CH3OH, and CO2 mixed in the 
ice. At high CH3OH concentrations, both CH30H:C02 
and H20:CH30H:C02 laboratory ices have a ^^C02 band 
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Fig. 6. a Observed peak position and 
FWHM of solid ^^cOs for all ob- 
served sources. The lines represent 
the peak and FWHM of the labo- 
ratory experiments as presented in 
Fig. |2[ The dashed line is the mix- 
ture H2O:CH3OH:CO2=l:0.92:l at 
the indicated temperatures, b Same 
as panel (a), but here the observed 
interstellar ^^C02 band is decom- 
posed with Gaussians. This Figure 
shows that most interstellar ^'^C02 
bands can be explained by a sin- 
gle mixture at different temperatures 
(a) or by a strongly heated "segre- 
gated" ice component together with 
a non-segregated (perhaps cold) com- 
ponent (b). 
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Fig. 5. Comparison of the spectra of Elias 16 and the K gi- 
ant a Bootis on optical depth scale. The bottom spectrum 
is the subtraction of a Bootis from Elias 16, to correct for 
photosphcric lines. The * symbols on top indicate the ex- 
pected frequencies of the gas phase CO band heads (see 
text). The residual feature at '^2278 cm~^ is attributed 
to ^^C02 ice in the Taurus molecular cloud. 

that peaks at frequencies less than ~2276.0 cm~^. Thus, 
we can exclude that the ice is CHaOH-rich, i.e. we con- 
clude that CH30H/C02< 1, which is in good agreement 
with the observed column densities (Sect. 4.4). 

For several reasons, the broad component of the ^^C02 
absorption band does not provide further constraints on 
the CH3OH/CO2 mixing ratio, nor the H2O/CO2 ratio. 



First, although the width of the absorption band is a good 
discriminator between polar and and apolar ices (Sect. 3), 
it depends only weakly on the relative abundance of H2O 
and CH3OH in the ice (Fig. ^). Second, for peak positions 
in the range 2277.5-2279.5 cm~^ there is an ambiguity be- 
tween H2O, and CH3OH ices. Good fits to the interstellar 
spectra can be obtained with H20:C02, CH30H:C02, as 
well as H20:CH30H:C02 ices at a range of temperatures. 
Thus, although we can put an upper limit of CH3OH/ 
CO2 < 1, no lower limit to the CH3OH concentration 
can be set from the broad ^'^C02 component. Further- 
more, for most objects good fits are obtained for a wide 
range of H2O concentrations. In some lines of sight we can 
put a significant lower limit H2 0/C02> 1, when fitting 
H20:C02 ices (e.g. Elias 29, NGC 7538 : IRS9, W 33A). 
However, when adding CH3OH to the ice, equally good 
fits with even lower H2O concentrations can be obtained. 

The peak position of the narrow component is remark- 
ably constant with a weighted mean of 2282. 5±0. 2 cm^^, 
while the width is typically FWHM=2.5 cm^^. This small 
width can be ascribed to ^'^C02 interacting with apolar 
molecules (Figs. ^ and ||b). The accurate peak position 
provides further constraints to the origin of this feature. 
For some sources (GL 2136, W 3 : IRS5), both pure CO2 
and heated "segregated" polar CO2 ices provide good fits. 
However, for the sources with the strongest narrow peaks 
(S 140 : IRSl, and GL 2591), the best fits are only obtained 
by heated polar CO2 ices, showing segregation behavior. 
These ices must contain CH3OH, and may contain H2O 
(but not necessarily). The initial CO2 concentration in 
these heated CHsOH-containing ices must be in between 
20-90%. At higher CO2 concentrations, the peak position 
and width are not well matched (Fig. At very low CO2 
concentrations, no segregation takes place upon heating 
(see Sect. 3 and the CH3OH:CO2=l:10 ice in Fig. |). 
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Table 2. The observed peak position (v), width (FWHM), and peak optical depth (r) of the interstellar solid ^■^C02 
band. For sources with high quality spectra we give in each second line the results of a decomposition of the profile, 
using 2 Gaussians. The standard deviation a is given in parentheses. 



Object 


u 


FWHM 


r 


V 


FWHM 


T 




















o -t Ar\ TT> O 1 

b 140 : IKbl 


2281.5 (0.8) 


5 (3) 


0.050 (0.005) 












2278.1 (1.8) 


9 (3) 


0.021 (0.003) 


2282.0 (0.6) 


3.0 (1.5) 


0.034 


(0.002) 


GL 2591 


2282.0 (0.8) 


7 (3) 


0.034 (0.005) 












2277.4 (1.5) 


8 (3) 


0.015 (0.002) 


2281.8 (0.3) 


2.6 (1.2) 


0.028 


(0.003) 


GL 2136 


2282.5 (0.8) 


10 (1) 


0.071 (0.008) 












2276.7 (0.5) 


7 (1) 


0.056 (0.004) 


2282.5 (0.5) 


2.5 (0.6) 


0.055 


(0.003) 


TXT" O TTlOtr 


2283.0 (0.8) 


8.9 (0.8) 


0.073 (0.012) 












2279.0 (0.8) 


8 (2) 


0.050 (0.008) 


22v>6.2 (0.3) 


1.3 (0.7) 


0.049 


(0.005) 


JNGG 7538 : msQ 


2279.0 (2.0) 


10.4 (0.8) 


rxir ( r\ r\-\ \ 

0.15 (0.01) 












2277.6 (0.2) 


8.8 (0.4) 


0.138 (0.004) 


2282.6 (0.2) 


3.0 (0.4) 


0.091 


(0.003) 


Elias 29" 


2279.0 (2.0) 


10 (1.5) 


0.071 (0.009) 












2277.8 (1.2) 


8.6 (2.3) 


0.064 (0.004) 


2282.3 (1.4) 


2.3 (2.9) 


0.026 


(0.005) 


W 33A" 


2276.5 (1.0) 


9.5 (0.8) 


0.225 (0.015) 












2275.3 (0.2) 


6.4 (0.4) 


0.177 (0.008) 


2280.0 (0.4) 


7.2 (0.7) 


0.119 


(0.008) 


Elias 16 


2278.0 (2.0) 


8.5 (1.5) 


0.060 (0.017) 










GC 3 


2278.0 (0.8) 


4(2) 


0.042 (0.015) 










R CrA : IRS2" 


2280.0 (1.0) 


7(1) 


0.13 (0.02) 












2278.4 (0.6) 


6.9 (1.1) 


0.092 (0.008) 


2281.3 (0.6) 


3.8 (1.2) 


0.064 


(0.005) 


NGC 2024 : IRS2 


2278 (2) 


8(2) 


0.060 (0.015) 










HH 100" 


2278 (2) 


8.5 (1.0) 


0.100 (0.015) 












2277.5 (0.4) 


5.9 (0.7) 


0.098 (0.008) 


2281.9 (0.5) 


3.1 (1.0) 


0.053 


(0.005) 


GL 989'' 


2281 (2) 


10 (3) 


0.046 (0.008) 











^Equally good fits are obtained with 1 and 2 Gaussians (see text). 
''Moderate spectral resolution; narrow component unresolved. 
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Fig. 7. Laboratory fits to the narrow CO2 component 
of S 140 : IRSl, indicating that a heated "segregated" 
CH3OH containing ice, rather than a pure CO2 ice is 
needed to fit the exact peak position. 



narrow features in GL 2591, and S 140 : IRSl are also 
obtained by CO-rich ices, peaking at 2281.7 cm^^. How- 
ever, in these l ines o f sight no CO ic e was detected (van 
Dishoeck et al. |l996[ Gerakines et al. |l999| ; Table |), and 
these mixtures can be excluded. 

Concluding, the detected narrow 2282 cm~^ compo- 
nent can be well fitted with a heated polar CO2 ice show- 
ing segregation. For at least two sources in our sample, 
this ice must contain CH3OH (CH3OH/CO2>0.10) to fit 
the exact peak position. For the other sources, a pure CO2 
gives reasonable fits as well, but mixtures with significant 
amounts of other apolar species (CO, O2) can be excluded. 



4.2.2. Single component ices 

When comparing the observed peak frequencies and over- 
all widths (not fitting Gaussians) with the laboratory re- 
sults, we find that toward most lines of sight the com- 
plete "'^'^002 ice band profile can be fitted with one single 
ice at a specific laboratory temperature (Figs. ||a and ||). 
Good fits are obtained with an ice in which H2O, CH3OH, 
and CO2 are about equally abundant, which is in good 
agreement with studies of the ^^C02 bending mode (Ger- 



For completeness, we mention that satisfactory fits to the akines et al. 1999). However, equally good fits are obtained 
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with other mixing ratios as well. For objects showing ev- 
idence for a separate 2282 cni^^ component (e.g. S 140, 
GL 2591), the ice must contain at least some CH3OH to 
fit the peak position (CH3OH/CO2>0.10) and the H2O 
abundance must not be too large (H20/C02<3). To fit 
W 33A with a single ice, the CH3OH abundance must be 
at least 30% of CO2. In contrast, the object R CrA : IRS2 
must have an CH3OH/CO2 ratio less than ~50% (note 
the offset from the other sources in Fig. ||a). For other 
objects, such as HH 100, the ice composition can not be 
significantly constrained. 

If we do assume that the "'^'^002 band profile origi- 
nates from a single ice, and it contains CII3OII and II2O 
with comparable concentrations, the temperatures of the 
best fitting laboratory spectra to the sources showing the 
2282 cm~^ substructure are well constrained within the 
range T = 115 — 118 K (Figs. ^, and ||; Gerakines et al. 
1999| ). At the long time scales in interstellar space, the 



segregation takes place at a much lower temperature of 
-77 K (Sect. 5.1). 

4-.3. CO2 column densities 

Column densities of solid "'^'^C02 are given in Table ^ and 
were derived by integrating the optical depth in the fre- 
quency range 2269-2288 cm~^ and dividing this by the 
integrated band strength A(^^C02). Laboratory experi- 
ments indicate that yl(^^C02) depends on the ice compo- 
sition, but only weakly on temperature (Gerakines et al. 
1995| ). The band strength is 7.8 10"^'' cm molecule"^ for 
pure CO2 ice, and similar for the mixture Il20:C02=1.6:l, 
but 15% lower for H20:C02=24:1. Since in the interstel- 
lar medium typically H20/C02< 5 (Sect. 4.4), we take 
yl(^^C02)=7.8 10"^^ cm molecule"^ Note that the varia- 
tion of A(^^C02) seen i n apolar CO and O2 containing ices 
(Gerakines et al. 1995) is not applicable, since our study 
shows that interstellar CO2 is absent in these ices. The er- 
ror bars were determined from the average signal-to-noise 
values given in Table |l]. 

Using these ^'^C02 column densities, we are for the 
first time able to calculate the ^^C02/^^C02 abundance 
ratio from the solid state, in a large number of sight-lines. 
The "'^^C02 column densities were derived from the ^^C02 
bending and stretching modes (~4.27 and 15.2 ^m; Ger- 
akines et al. 1999). The column densities derived from 



both modes agree very well within the error bars given 
in Gerakines et al. (|l999|) . These errors include ISO-SWS 
calibration uncertainties (mainly applicable to the stretch- 
ing mode) and uncertainties due to the continuum de- 
termination, which is particularly difhcult for the ^^C02 
bending mode, since it is blended with the bending mode 
of silicates (Gerakines et al. 199£; Boogert 1999 ). 

The ^^C02/^^C02 column density ratio has a rela- 
tive error less than 16% in 6 lines of sight (Table 
However, laboratory experiments have shown that the 
^^C02/"'^'^C02 ratio can not be determined with an ac- 
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(cm 
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Fig. 8. Laboratory fits to the interstellar ^^C02 band. 
For each source, the left panel gives a two component fit 
with the thin solid line representing a mixture of CO2 
with II2O or CII3OII. The dashed line is the mixture 
H20:CH30H:C02=1:1:1, heated to T =125 K. Note the 
significant variations in the relative contribution of these 
components. The thick gray line is the sum of both labo- 
ratory spectra. The right panel shows for each source the 
best fit with H2O:CH3OH:CO2=l:0.92:l at the indicated 
temperature. For sources that do not give good fits (e.g. 
HH 100), changing the mixing ratio somewhat will im- 
prove the fit at the indicated temperature. The luminous 
protostars in the left panels are ordered with a decreasing 
strength of the narrow 2282 cm^^ peak. 
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Table 3. The CO2 column density and the 
ratio derived from solid CO2 and gas phase CO. 



12 



Object 


iV(13 


CO2) 




13 c 






lO^'^cm-^ 


CO2 


CO'' 


kpc 


GC 3 


0.21 


(0.04) 


52 (11) 


24 (lY 


0.5 


W 33A 


2.74 


(0.21) 


53 (8) 


39 (I)'' 


4.5 


GL 2136 


0.73 


(0.05) 


107 (8) 




6.1 


Elias 29 


0.83 


(0.05) 


81 (11) 




7.8 


R CrA : IRS 2 


1.17 


(0.09) 


73 (16) 




7.9 


HH 100 


1.03 


(0.08) 


52 (11) 




7.9 


GL 2591 


0.26 


(0.03) 


62 (10) 




7.9 


Elias 16" 


0.47 


(0.15) 


98 (38) 




8.2 


S 140 : IRSl 


0.38 


(0.03) 


111 (9) 




8.4 


NGC 2024 : IRS2 


0.55 


(0.13) 


105 (33) 




8.4 


GL 989 


0.62 


(0.09) 


84 (21) 


62 (3) 


8.7 


NGC 7538 : IRS9 


2.03 


(0.12) 


80 (11) 




9.4 


W 3 : IRS5 


0.63 


(0.13) 


113 (37) 


66 (4)" 


9.7 



^ ^''COg band corrected for photospheric OH 
absorption lines 



taken from Langer & Penzias (1990) 
" Sgr B2, at (ARA, ADec)=(14.7', 26.1 

W 33(0,0), at (ARA, ADec)=(-6.0', 
« W 



at (ARA, ADec)= 



from GC 3 
3.8') from W 33A 
(10.5', -13.3') from 



3(0H) 
W 3 : IRS5 

* galacto-centric radius calculated from Galactic coordinates 
from Simba d database and distances mentioned in Aitken 
19931 (GL 989, GL 2136, W 3 : IR S5), Chernin [ggf 
IRS2), Giirtler et al. |l99l| (GL 2591, S 140 



et al. 



et al. 1990 (GC 3) 



L996 



(NGC 2024 

NGC 7538 : IRS 9, W 33A), Wliittet~|l974| (Elias 29), Sara- 
ceno et al 



(HH 100, R CrA : IRS2, Elias 16), Okuda 



curacy better than 20% (Sandford & AUamandola 1990; 
Ehrenfreund et al. 1997), unless the ice composition is ac- 
curately known. The ^^C02/^^C02 ratio of integrated 
optical depths is 19% lower in an H2O:CO2=l:10 ice 
compared to pure CO2, and 14% higher for H2O: CO2— 
5:1. For 3 objects with low statistical errors (GL 2136, 
GL 2591, and S 140), we find that the CO2 is heated 
and segregated, closely resembling the band profile of pure 
C02- Thus, the band strengths are likely similar to pure 
CO2. However, for Elias 29, W 33A, and NGC 7538 : IRS9 
we can not exclude that a large fraction of the CO2 is em- 
bedded in an H20-rich ice, and we conservatively assume 
a 14% error in the ^^C02/^^C02 column density ratio for 
these sources. 

Interesting variations in the ^^C02/^'^C02 ratio are 
evident. Most objects have ratios in between 70-110, but 
some are significantly lower, e.g. W 33A, GC 3 and 
HH 100. We will discuss these results in Sect. 5.3. 



C/ C 4-4- Further constraints on the solid H2O/CO2 and 
CH3OH/CO2 ratios 

The composition of interstellar CO2 ices can be further 
constrained by the CH3OH, H2O, and CO2 column den- 
sities along the observed lines of sight. The H2O/CO2 ice 
column density ratio is typically ~ 5, but can be as high as 
8 (G L 2136 , GL 2591) and as low as 1.5 (GC 3; Gerakines 
et al. [19991) . Thus, although mixtures with H2O/CO2>10 
do provide good fits to most of the observed sources (e.g. 
W 33A, R CrA : IRS2), these mixtures are unrealistic. 
Furthermore, observations of the C-H stretching mode of 
CH3OH (3.54 /^ m) in dicate that CH3O H/H2O<0.10 (Al- 
lamandola et al. [1992 ; Chiar et al. 1996| ), resulting in typi- 
cally CH3OH/CO2<0.50. For some individual sources, the 
observed CH3OH/CO2 ratio is higher (W 33A: CH3OH/ 
C02= 1.90), or lower (Elias 16: CH3OH/CO2 <0.20). 
This is consistent with our conclusion that the interstellar 
^^C02 ice band for most sources is best fitted with labora- 
tory ices with a low CH3OH concentration, i.e. for which 
O.K CH3OH/CO2 < 1. To further limit the CH3OH/CO2 
ratio in the ice, it is, in particular for the low mass pro- 
tostars, important to obtain very sensitive observations of 
the C-H stretching mode of CH3OH. 



5. Discussion 

5.1. Heating of CO 2 ices 

The ^^C02 ice band profile varies remarkably between the 
observed lines of sight. Most significant is the detection of 
a narrow substructure at ^2282 cm^^ toward the mas- 
sive protostellar objects S 140 : IRSl, GL 2136, GL 2591, 
and W 3 : IRS5. None of the low-mass protostars show 
this component, nor do some high-mass objects (W 33A, 
NGC 2024 : IRS2), the Galactic Center object GC 3, and 
quiescent cloud material traced by Elias 16. Laboratory 
simulations indicate that this distinct profile can be at- 
tributed to a heated, "segregated" polar CO2 ice. Here, 
we further investigate this hypothesis by comparing the 
strength of the narrow "'^'^C02 peak, normalized to the 
Gaussian depth of the broad component, to a number of 
established temperature tracers (Table ||). 

Although only sources with a bolometric luminosity 
Lboi > Lq have detected narrow ^^C02 peaks (Ta- 
ble ^; Fig. H), a high luminosity can not be the only nec- 
essary condition. There is no, or very weak, sign of the 
^^C02 substructure in the luminous objects W 33 A and 
NGC 7538 : IRS9. The dust temperature in these sight 
lines may be too low to cause the ice to segregate. The 
interstellar CO2 ice afi^ected by segregation must have a 
temperature of T = 50-90 K (or ~100-180 K in the labo- 
ratory), depending on whether the ice is H2O, or CH3OH- 
rich. Assuming the dust radiates as a blackbody, modi- 
fied by a power law emissivity with index —1, this corre- 
sponds to emission peaking at wavelengths of ~25-50 fj,m. 
Non-heated ice, at temperatures of ^20 K, radiates near 
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Fig. 9. a— c The strength of the narrow ^^C02 peak plotted against various physical parameters of the protostars. High 
and low luminosity objects are indicated by bullets and triangles respectively. Panel a plots the bolometric luminosity, 
showing a poor correlation. The ratio of hot over cold dust emission in panel b gives a much better correlation. The 
far infrared flux ratio for W 3 : IRS5, indicated with an open symbol, is unreliable due to contamination in the large 
beam. Panel c shows that there is also a fairly good anti-correlation with the solid CO2 abundance. 



100 /im. Hence, a useful parameter for comparison is the 
flux ratio 



hot 



F(45 /m) 
" i^(100 ^m) 



(1) 



where the F{A5 /im), and F(100 /^m) fluxes are deter- 
mined from ISO-Long Wavelength Spectrometer (ISO- 
LWS) spectra. We checked for contamination by extended 
emission in the 80" ISO-LWS beam, by comparing with 
the ISO-SWS flux at 45 nm in a 20" beam. For most 
sources there is good agreement within the instruments 
cross calibration uncertainty of ~ 30%. An exception is 
W 3 : IRS5, with a difference of more than a factor of 
2, and we consider the flux ratio unreliable. We find that 
sources with a large narrow ^^C02 peak strength have a 
high hot/cold dust flux ratio (Table |; Fig. ^. And indeed, 
sources with low upper limits to the ^■^C02 substructure, 
such as W 33A, NGC 7538 : IRS9, and Ehas 29, are sur- 
rounded by much larger fractions of cold dust. 

A similar picture is seen in the gas phase characteris- 
tics. Near-infrared CO observations indicate that a large 
fraction (~50%) of the gas toward the objects with the 



narrow 



i3r 



^C02 ice substructure is hot, T = 400 - 1000 K 
(Table It must be noted that in all objects the temper- 
ature of the hot CO gas is well above the ice sublimation 
temperature T = 90 K, and thus the observed heated CO2 
ice originates from a region outside the hot core. On the 
other hand, the temperature of the cold CO gas compo- 
nent is in some lines of sight too low to cause the CO2 
ice to segre gate (e.g. T = 28 K for S 140 : IRSl; MitcheU 
et al. 1990). Perhaps the observed hot CO2 ice originates 



from the interface between the hot core and the cool sur- 
rounding medium. Indeed, other gas phase temperature 
tracers (e.g. CH3CN, CH3OH; van der Tak et al. [1999D 



reveal temperatures that are in between the cold and hot 
components found by CO studies. 

Finally, the abundance of volatiles, such as ices, is ex- 
pected to be a good tracer of the dust temperature along 
the line of sight. Indeed, sources with low CO2 abun- 
dances have the strongest substructure in the ^^^002 ice 
band (Table ^; Fig. ^). For solid CO, a similar trend is 
observed. The objects with the strongest narrow ^'^C02 
peaks, S 140 : IRSl and GL 2591, have no detected CO 
ice along the line of sight. This contrasts with W 33 A and 
NGC 7538 : IRS9, which even have detections of highly 
volatile apolar CO detections (Smith et al. 1989; Tielens 
et al. 1991). All low-mass objects have a high fraction of 



apolar CO, and at the same time show no sign of heated 
CO2 ice along the line of sight. 

Thus, we flnd that the strength of the narrow '^'^002 
component at 2282 cm~^ correlates with the dust and CO 
gas temperature along the line of sight, and anti-correlates 
with the ice abundances. This is good evidence that the 
observed feature can be attributed to a heated polar CO2 
ice, and it confirms our results from the laboratory fits. 
Thermal processing is thus an important process in in- 
terstellar space, at least toward luminous protostars. The 
differences in the various lines of sight are intriguing, and 
may indicate an evolutionary sequence with increasing 
thermal processing. Perhaps several distinct layers exist 
around massive protostars, in a way sketched in Ehren- 
freund et al. (1998): a hot core where the ices have evap- 
orated, a warm region surrounding it, where ice crystal- 
lizes, and further out a colder region. With time, these 
temperature regions expand outward. This would corre- 
spond to the two component ice model, where the narrow 
peak and broad components arise in physically different 
regions (Sect. 4.2.1). 
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Table 4. Strength of the narrow ^■^C02 ice peak and observed physical parameters from hterature 



Object 


narrow 13*^/^ ^ 
broad '-''-'2 


-^hot CO 


Thot CO 




CO 


Lhoi 










% 


K 


% 


ref. 




ref. 




-t Ar\ TT> 1 

b 140 : IKbl 


1.6 (0.3) 


41 (10) 


390 





[4] 


3 10* 


[10] 


3.2 


GL 2591 


1.9 (0.4) 


56 (10) 


1010 





[2] 


6 10* 


[10] 


4.4 


GL 2136 


1.0 (0.1) 


32 (16) 


580 


< 10 


[3] 


7.2 10* 


[9] 


2.8 


W 3 : moo 


1.0 (0.1) 


50 (25) 


580 


80 




4.7 10^ 


[15] 


/I 1 / 
4.1 


NGC 7538 : IRS9 


< 0.7 


2 (1.5) 


180 


92 (11) 


[1] 


9.2 10* 


[14] 


2.1 


Elias 29 


< 0.5 






88 


[5] 


36 


[13] 


2.0 


W 33A 


< 0.5 


52 (26) 


120 


26 (3) 


[1] 


1.1 10^ 


[10] 


1.3 


Elias 16 


< 0.4 






86 (4) 


[6] 








GC 3 


< 0.3 









[7] 








R CrA : IRS 2 


< 0.7 






92 (34) 


[1] 


12 


[8] 




NGC 2024 : IRS2 


< 0.9 






74 


[3] 


1 10^ 


[12] 




HH 100 


< 0.6 






63 (7) 


[1] 


14 


[8] 


1.1 


GL 989 


e 






71 


[3] 


3.3 10^ 


[11] 





^ narrow over broad Gaussian ^^C02 peak depth ratio 

^ hot CO gas column density in percentage of total N{CO gas); Mitchell et al. 1990| 

column density of CO in apolar ice in percentage of total A'^(CO ice). Entries with '0' indicate 

that no CO ice was detected, with a significant upper limit. 

calculated as F(45^m)/_F(100/im) from ISO-LWS spectra 

° line profile unresolved 

f 



observed FIR flux heavily contaminated by nearby objects 
1] Chiaret al 
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Alternatively, the very good fits that we obtain with 
single ices, would suggest that all the ice has one spe- 
cific temperature, and thus is present in a specific re- 
gion around the protostar. The temperatures of the least 
(W 33A; T=115 K) and most (S 140; r=118 K) evolved 
luminous objects are extremely well constrained in an in- 
terval of only 3 K. At the long time scales in interstel- 
lar space, this small temperature difference can be trans- 
lated to a difference in heating time at the amorphous- 
crystalline phase transformation (segregation) tempera- 
ture. Amorphous ice has a highly disordered structure, 
characterized by a broad distribution of bond angles be- 



tween neighboring H2O molecules (Madden et al. 1978). In 
contrast, crystalline ice has well defined bond angles. This 
broad distribution implies a range of activation energies 
associated with the amorphous-crystalline phase transi- 
tion. At a relatively low temperature in the laboratory 
(r = 115 K), only the lowest barriers can be surmounted. 
As the temperature increases (to T = 118 K), the higher 
barriers can be scaled as well. Hence, this temperature 
range implies a variation in barrier height by about 3%, 
with a typical height of i?sGgr=4600 K (calculated from a 
5500 K ba rrier at 140 K for pure ice; Tielens & AUaman- 
dola 1987). In interstellar space, rather than temperature. 



time is the essential parameter. The time scale Tsogr needed 



to surmount an energy barrier i^sogr at temperature T is 
given by: 



' scgr — ^0 



(2) 



with 



1^0 



5 10 s the 0-H ben ding m ode 



bration frequency (Tielens & AUamandola 1987 ). As the 
protostar is formed, it heats the surrounding molecular 
cloud, creating a hot core region around it where the ice 
evaporates. Immediately surrounding this region is a zone 
where the ice has been warmed sufficiently to start the 
amorphous- crystalline phase transformation, correspond- 
ing to the r = 115 K case in the laboratory. As time 
progresses, this phase transformation progresses as well, 
and in essence, although the temperature does not change, 
the interstellar ice corresponds to warmer and warmer (i.e. 
115-118 K) laboratory experiments. At the typical time 
scale of evolution for hot core regions of r ^3 10^ years 



(Charnley et al. 1992), the laboratory temperature range 
of 115-118 K effectively corresponds to a segregation tem- 
perature of ~77 K (Eq. 2). Thus, at this temperature, 
W 33A and NGC 7538 : IRS9 would have the youngest 
and least evolved hot cores, while S 140 and GL 2591 are 
~ 3 10^ years older and have succeeded to heat the sur- 
rounding gas and ice significantly. However, the derivation 
of a heating time scale is complicated by the presence of 
temperature gradients, and possible intrinsic differences 
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between the sources. It seems thus hkely that both time 
and temperature effects play a role in the evolution of in- 
terstellar CO2 ices. Detailed modeling of the temperature 
and density structure around massive protostars is needed 
to derive evolutionary time scales (e.g. Van der Tak et al. 



199c). 



5.2. Chemistry of interstellar CO 2 

Interstellar CO2 may form on grain surfaces, through the 
reaction of atomic O with CO, atomic C, or C+, accreted 
from the gas phase. Although CO, rather than C or C"*", is 
generally assumed to be the origin species, solid and gas 
phase isotope ^^C/^^C ratios provide a good test for this 
(Sect. 5.3). In this grain surface scenario, the ice compo- 
sition provides a measure for the gas phase abundances 
during accretion, convolved with reaction efficiencies. The 
ratio of CII3OII/CO2 is a measure for hydrogenation over 
oxidation reactions of CO. The ^'^C02 profile provides an 
upper Hmit CII30II/C02<1 in the ice, and thus oxidation 
of CO is at least as efficient as hydrogenation, despite the 
much larger atomic H abundance compared to atomic O. 
At the same time, the solid H2O/CO2 ratio measures the 
gas phase (0-|-02)/C0 ratio at accretion, since II2O is 
formed from O2 or O3. The ^'^C02 profile does not sig- 
nificantly constrain the solid H2O/CO2 ratio, it could be 
as high as the ra tio of column densities, i.e. H20/C02<5 
(Gerakines et al. 1999 ). Further constraints to the chemi- 
cal conditions during CO2 formation on interstellar grains 
can be obtained from th e ^^C 02 bendi ng an d stretching 
modes (Gerakines et al. 1999 ; Boogert 1999 ). In partic- 
ular, the prominent red wing in the bending mode is a 
signature for the presence of C02:CIl30II clusters in in- 
terstellar ices. 

An alternative way to form interstellar CO2 is through 
irradiation of CO containing ices by ultraviolet (UV) pho- 
tons. Even in "cold" sight-lines such as Elias 16, the UV 
flux from the ISRF or induced by cosmic rays may be 
high enough to explain the observed CO2 abundance (see 
Whittet et al. 1998 for a quantitative discussion). The oxy- 
gen atoms are liberated from H2O in polar ices, or O2 in 
apolar ices. The interstellar ^^C02 profile does not show 
the narrow signatures of CO2 mixed with apolar molecules 
such as CO and O2. Hence, if UV processing is the dom- 
inant mechanism, CO2 most likely originates from polar 
CO ices. Indeed, interstellar CO is presen t in s ignificant 
quantities in the polar phase (Tielens et al. 1991 Chiar et 
al. 119981) . 



5.3. Galactic ^^C/^^C abundance ratio 

The derived ^^C02, and published ^^C02, column den- 
sities were used to calculate ^^C/^^G abundance ratios 
(Table |). It is the first time that the ^'^C/'^^C ratio is 
determined from the solid state. Previously, gas phase 
molecules were used for this purpose, using UV, optical. 



and radio observations (see Wilson & Rood 1994 for an 
overview). C^^O and H2CO observations yield an increas- 
ing ^^C/^^C ratio with Galacto-centric radius (Langer & 
Penzias |1990D . This can be understood by the different ori- 
gin of ^^C and ^^C atoms. The ^^C is rapidly produced 
by Helium burning in massive stars, and injected into the 
interstellar medium by supernovae. In later stellar gener- 
ations, ^^C seeds are converted to ^^C in the CNO cycle 
of low- and intermediate-mass stars, during their red giant 
phase. This is a much slower process. Thus, the ^^C/^'^C 
ratio is a measure for the enrichment of the interstellar 
medium by primary to secondary stellar processes. The 
observed gradient in the Galactic ^^C/^'^C ratio can then 
be understood by a higher star formation rate in the inner 
parts of the Galaxy (Tosi |l982[ ). 

The ^^C/^^C ratio derived from solid CO2 indicate 
the same trend as found in gas phase studies (Table |[ 
Fig. no). Low values are found near the Galactic Cen- 
ter (^02/"C02=52±ll; GC 3) and in the molecu- 
lar ring, at ~4 kpc from the Galactic Center (W 33A; 
53±8). Objects at larger Galactic radii mostly have higher 
^^002/^^002 ratios. We have fitted a linear relation to 
the ^^002/^^002 ratio as a function of Galacto-centric 
radius, where, as in the gas phase studies, we exclude the 
Galactic Center from the fit. We find the following rela- 
tion, using the errors as statistical weights: 



(12C02/^^C02) = (4.5 ± 2.2)i?G + (48 ± 16) 



(3) 



If instead we take an unweighted fit, like in some gas 
phase studies (Wilson & Rood 1994), we find a gradi- 
ent of 6.5±4.5 and absisca 33±21. Thus, the evidence for 
a gradient in the ^^C/^'^C ratio derived from solid CO2 
is weak. Within the (rather large) error bars, this gradi- 
ent is comparable to the H2CO (6.6±2.0) and CO data 
(7.5±1.2). It must be mentioned that a limitation of our 
analysis is the lack of solid CO2 observations at low Galac- 
tic radii (3-6 kp c). W e note that also in gas phase studies 
(Wilson & Rood 1994), there is considerable scatter in the 
12q^13(^ Yaiio, and no clear trend at Galactic radii > 6 kpc 
is present either (Fig. p^ ). 

Although the gradients agree well for CO2, CO, and 
H2CO, we find that ^^C/^^C ratios derived from CO (ab- 
scissa — 2.9±7.5 in Eq. 2) are systematically somewhat 
smaller. A comparison of the gas and solid state ^^C/^'^C 
ratios for individual objects is possible for only 4 sources. 
For all of these, our solid CO2 observations yield higher 
12q^13(^ ratios compared to gas phase CO observations 
(Table H). However, for both the overall linear fit (Eq. 2) 
and for these individual sources, the differences are statis- 
tically small (< 3(7). For a proper interpretation of these 
possible differences it is important to note that the gas 
phase observations were done several arcminutes away 
from our CO2 observations (Table In the gas phase, 
considerable variations of the ^^C/^^C ratio were reported 
within the same region, such as within the W 33 cloud 
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2468 10 0246 
Galacto-centric Radius (kpc) 

Fig. 10. The ratio of solid ^^002 and ^'^C02 column den- 
sities, versus the Galacto-centric radius of the observed 
sources (left panel). The solid line is a linear weighted fit 
to the data. The right panel shows the ^^C/^'^C ratio from 
gas phase studies (bullets: C^^O, open circles: H2CO) as 
taken from Wilson & Rood (1994). The dotted and dashed 
lines are linear fits to these gas data. 



(Langer & Penzias 1990D and toward the supernova rem- 
nant Gas A (Wilson & Rood 1994). In particular, compar- 
ison of our ^^G02/^'^G02 ratio toward the Galactic Genter 
object GG 3 with the gas phase determination toward the 
Galactic Genter supernova remnant Sgr B2 is complicated 
by the quite different stellar evolution history that these 
regions may have. 

Finally, we determined that i2c02/i3G02=69±15 
(1 a) for the local ISM, by taking an unweighted mean for 
the three objects located in the p Ophiuchi, and Gorona 
Australis clouds. We excluded Elias 16 from this mean, 
since it has a much larger error bar than the other sources 
(Table ||). This value is comparable to the ^^G/^'^G ratio 
in the local ISM, as determined from atomic G (58±12; 
Keene et al. |1998[ ) and G+ (58±6; Boreiko & Betz |l996D 
observations, as well as from G^^O and H 2GO observations 
(12c/ i3g=77±7; Wilson & Rood |l994[ ). 

In summary, we find that the gradient of the ^^G/^^G 
ratio with Galacto-centric radius is the same for solid 
GO2 and gaseous GO and II2GO. In this relation, the val- 
ues derived from GO, however, tend to be smaller com- 
pared to GO2 and H2GO. The value for the local ISM 
(^^G02/^'^G02= 69±15) is in reasonable agreement with 
gas phase studies. These results could have important im- 
plications for our knowledge of interstellar chemistry. Both 
gaseous GO and H2GO may have been affected by chem- 
ical fractionation and isotope selective destruction (e.g. 
Langer et al. 1984; Tielens 1997). If fractionation is an 
important mechanism, preferential incorporation of -'^^G"*" 
in GO, leads to low i^CO/^^cO and high ^^C+ /^^C+ ra- 



tios. On the other hand, if isotope selective destruction 
dominates over fractionation, ^■^GO is preferentially de- 
stroyed, leading to high i^CO/^^GO and low ^^C/^^C ra- 
tios. Garbonaceous molecules derived from G"^ or G would 
then have high or low ^^G/^^G ratios, depending on the 
dominant process and chemical pathway. Our lines of sight 
trace dense molecular clouds, and GO is the main reser- 
voir of G. Then, fractionation or selective destruction will 
highly infiuence the i2c(+)/i3c(-i-) j-^^io, but will hardly 
affect the i^CO/^^GO ratio. The generally higher H^^cQ/ 
Hj^GO ratios compared to G O thu s suggest that H2GO 
originates from 0+ (Tielens |l997 ). Similarly, if indeed 
the ^■^002/^^002 ratios are higher compared to GO, this 
would imply that GO2 is formed from G"*" rather than 
G or, as is generally assumed, from GO (van Dishoeck 
et al. 1996 ). However, we must emphasize that with the 
present large uncertainties in the ^^G/^^G ratios we can 
not make strong statements about the chemical pathway 
to form interstellar GO2. We expect that with future im- 
provements of the ISO-SWS data reduction, the large er- 
ror bars in the GO2 abundances for some sources (GL 989, 
NGG 2024 : IRS2, and W 3 : IRS5) can be reduced. Fur- 
thermore, it is important to obtain ^^GO/^'^GO ratios in 
the same line of sight as our GO2 observations, preferably 
by absorption line studies along a pencil beam. Observa- 
tions of the solid ^^GO/^'^GO ratio would be particularly 
useful. This will also shed light on the origin of the large 
scatter seen in the ^^C/^^C ratio as a function of Galacto- 
centric radius for both solid GO2 (i.e. note the particularly 
high value for GL 2136 in Fig. 10) and gas phase species. 



It has been suggested that perhaps most of this scatter is 
real, and can be attributed to local variations in the star 
formation history (Wilson & Rood 1994| ). 



6. Summary and Conclusions 

We have presented ISO-SWS observations of the stretch- 
ing mode of ^'^G02 ice in the spectral range 2255-2300 
cm~^ (4.34-4.43 fim) in 13 Galactic lines of sight. All 
sight-lines show an absorption feature with a peak position 
in the range 2276-2280 cm-\ and a Gaussian FWHM of 
typically 6-9 cm~^. Additionally, the four high-mass pro- 
tostars GL 2136, S 140 : IRSl, GL 2591, and W 3 : IRS5, 
show a narrow (FWHM--3 cm-i) absorption fine at ~ 
2282.3 cm"-'^. This feature is much weaker or absent to- 
ward other high-mass objects (W 33A, NGG 7538 : IRS9), 
toward low-mass protostars, the Galactic Genter (GG 3), 
and quiescent molecular cloud material (Elias 16). 

These observational results are compared to labora- 
tory experiments of the stretching mode of ^^002 ice. We 
conclude that this band is a sensitive probe of the heat- 
ing history and composition of interstellar ice mantles. 
The profile of the interstellar band can be fitted in two 
different ways. First, the detected broad and narrow com- 
ponents could originate from ices with different heating 
histories and perhaps composition along the line of sight. 
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In this scenario, the detected broad component is ascribed 
to a non-hcatcd mixture of CO2 with polar molecules such 
as most likely H2O, and CH3OH. The mixing ratio of 
CH3OH/CO2 in the ice can generally be constrained to 
an upper limit of CH30H/C02<1, in agreement with the 
observed column densities along the lines of sight. The nar- 
row 2282 cm~^ absorption feature detected toward several 
high-mass protostars is ascribed to ^'^002 in a polar CO2 
ice that has been heated to temperatures of at least 50 K 
(in interstellar space). In such an ice, the bondings of the 
CO2, and polar molecules are segregated and the CO2 
band profiles resemble those in a pure CO2 ice. To fit the 
exact peak position of the 2282 cm~^ feature, the CO2 ice 
must contain at least some CH3OH (CH30H/C02> 0.1). 

A second way to fit the ^'^C02 band, is with a sin- 
gle ice at a specific temperature. This would locate the 
ice in a well confined region around the protostar. rather 
than in regions of different temperature as for the two 
component model. For the luminous protostars showing 
evidence for this narrow component, the best fitting labo- 
ratory ices have temperatures within a very small interval 
of T = 115 - 118 K (using H20:CH30H:C02=1:1:1). At 
the long time scales in interstellar space, this laboratory 
temperature interval corresponds to a heating time dif- 
ference between the luminous protostars comparable to 
the lifetime of hot cores (^3 10'* years) at an amorphous- 
crystalline phase transition temperature of ~77 K. How- 
ever, the derivation of heating time scales likely is com- 
plicated by the presence of temperature gradients, and 
possible intrinsic differences between the sources, such as 
different temperature and density structure. 

To further test the hypothesis that interstellar CO2 
ice is affected by heating, we calculate the ratio, and up- 
per limits, of the narrow to broad ^^C02 component peak 
depths. This quantity is compared to known physical pa- 
rameters of all objects. We conclude that the strength of 
the narrow ^^C02 substructure correlates with the dust 
and CO gas temperature along the line of sight, and anti- 
correlates with the ice abundances. This is further good 
evidence that the structure of ices around luminous proto- 
stars is affected by thermal processing, as concluded from 
the laboratory fits. This effect appears to be not impor- 
tant for low mass protostars and the quiescent medium. 
Although we have to keep in mind that our selection of 
low mass objects is small, this might imply that unaltered 
interstellar ices are included in comets. 

Finally, we have determined solid *''C02 column den- 
sities. These were used to derive ^^C/^^C ratios in all 13 
sight-lines. It is the first time that the Galactic ^^C/^^C 
ratio is determined from the solid state. Solid state de- 
terminations of the *^C/*'^C ratio are more reliable than 
gas phase studies, because the column density can be 
determined without uncertain radiative transfer effects. 
Like in gas phase studies, we find low values for the 
Galactic Center (GC 3; ^'^ CO ■2/'^^ CO 2 = 52±11) and the 
Galactic molecular ring (W 33 A; 53±8), and higher val- 



ues at larger Galacto-centric radii. The gradient of the 
^^C02/"'^^C02 ratio as a fmiction of Galacto-centric ra- 
dius (4.5±2.2 kpc~*) agrees with gas phase studies. The 
scatter in this relation is however large, and may perhaps 
indicate local differences in star formation history, since 
^^C is converted to *'^C in low and intermediate mass 
stars. Comparison of solid and gas phase isotope ratios 
can be used to trace the chemical history of molecules. In 
general, the *^C02/*'^C02 ratio, as for H2CO, tends to be 
larger compared to CO studies (~ 2.5 cr difference). In the 
dense molecular clouds that we have observed, where CO 
is likely unaffected by fractionation or isotope-selective de- 
struction, this would imply that interstellar CO2 is formed 
from C"'', rather than CO or C. However, to draw definite 
conclusions on interstellar chemistry and on global and 
local variations in the Galactic star formation history, we 
emphasize that a systematic determination of the *^C/*'^C 
ratio from atomic C^^\ CO, and solid CO2 in individual 
lines of sight is needed. At present, there is little over- 
lap between our solid CO2 observations and published gas 
phase observations. 
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